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We present for the first time a detailed and comprehensive analysis of the experimental results 
that set the current world sensitivity limit on the magnitude of the electric dipole moment (EDM) 
of the neutron. We have extended and enhanced our earlier analysis to include recent developments 
in the understanding of the effects of gravity in depolarizing ultracold neutrons (UCN); an improved 
calculation of the spectrum of the neutrons; and conservative estimates of other possible systematic 
errors, which are also shown to be consistent with more recent measurements undertaken with the 
apparatus. We obtain a net result of dn = —0.21 ± 1.82 x 10“^® ecm, which may be interpreted 
as a slightly revised upper limit on the magnitude of the EDM of 3.0 x 10“^® ecm (90% CL) or 
3.6 X 10“^® ecm (95% CL). 

This paper is dedicated by the remaining authors to the memory of Prof. J. Michael Pendlebury. 

PACS numbers: 13.40.Em, 07.55.Ge, 11.30.Er, 14.20.Dh 


I. INTRODUCTION 

Measurement of the electric dipole moment (EDM) of 
the neutron (or of other fundamental particles) provides 
an extremely sensitive approach to investigating poten¬ 
tial new sources of CP violation and physics beyond the 
Standard Model. 

The experimental technique underlying the most sensi¬ 
tive measurement to date mm of the neutron EDM has 
been discussed extensively in an earlier publication [3]. 
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The data were collected at ILL, Grenoble, between 1998 
and 2002. In this article, we focus solely on the analysis, 
which we have carried out anew. We begin, in Section 
[TI| with a description of the determination of the neutron 
resonant frequency via the Ramsey method of separated 
oscillatory fields in conjunction with the cycle-by-cycle 
corrections from the mercury co-magnetometer. 

A number of cuts were applied to the data. These are 
discussed in Section m 

Eor each data-taking run, of typically 1-2 days’ dura¬ 
tion, a value of the EDM was determined from the slope 
of a linear fit of the (field-drift corrected) neutron fre¬ 
quency versus applied electric field E. However, in the 
presence of a magnetic-field gradient, the mercury (and, 
to a far lesser extent, neutron) frequency acquires a com¬ 
ponent that is linear in the applied electric field [4], and 
thus mimics the signal for an EDM. We have no direct 
measurement of the applied H-field gradient, but since 
the neutrons have a lower center of mass than the (ther¬ 
mal) mercury atoms, any change in the gradient results 
in a change in the neutron-to-mercury frequency ratio. A 
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linear dependence of the measured apparent EDM as a 
function of this ratio therefore emerges, as shown in Fig. 
2 of [1] and reproduced here in Fig. with a positive 
(negative) slope for a downwards (upwards) direction of 
the applied magnetic field B. We discuss this further in 
Section nvi below. 




FIG. 1: (Color online) From [T]. Measured EDM (binned 
data) as a function of the relative frequency shift of neutrons 
and Hg. The solid red line is a linear best fit. 


The offsets and profiles of these lines can be affected by 
a number of factors, including nonuniform and/or hori¬ 
zontal magnetic-field gradients, the gravitationally en¬ 
hanced spectrum-dependent depolarization of the ultra¬ 
cold neutrons (UCN), and even the rotation of the Earth 
[5]. These are discussed in Section [V| 

In order to measure and compensate for such effects, 
measurements were made in an auxiliary bottle that used 
a lid of adjustable height, so that the neutron and mer¬ 
cury frequencies could be determined as a function of 
height for a given 5-field configuration. These measure¬ 
ments are described in Section ED 

In Section 1 ^ we describe how all of these measure¬ 
ments are brought together in a global fit, in order to 


determine the necessary corrections to the initial EDM 
estimate obtained from the fit to the two crossing lines. 

In Section VIII we discuss further possible contribu¬ 
tions to systematic errors. 

We close, in Section |IX[ with a final summary and 
conclusion. 

For convenience in what follows, and following the con¬ 
vention of [T], we define 


R 


a 


7Hg 

7n 


( 1 ) 


where z^, 7 are the respective frequencies and gyromag- 
netic ratios of the two species (neutrons and Hg), and we 
also introduce 


R' = Ra- 1. (2) 

In this context, R' is usually specified in parts per 
million (ppm). Throughout this analysis we have 
used values of 7 Hg = 7.590118(13) Hz//iT and 7 n = 
—29.1646943(69) Hz//iT, based on the measurements of 
Cagnac [ 6 ] and of Greene et al. [7]. Both of these mea¬ 
surements were made relative to the shielded proton gy- 
romagnetic ratio in pure water, for which the most re¬ 
cently measured value is 7 p = 42.5763866(10) Hz//iT 
[ 8 ]. From these results, we observe that R' = 0 when 
^n/^Hg = 3.8424560(66), which represents an accuracy 
of 1.73 ppm. We incorporate the recent measurement of 
7n/7Hg = —3.8424574(30) by Afach et al. [9] into the 
later stages of our analysis. 


II. NEUTRON FREQUENCY 
DETERMINATION 


The neutron frequency fitting procedure is described 
in m and [3]. For any measurement cycle (equivalently 
referred to as a “batch” in 121 ), the first-order estimate 
^n,Hg of tho neutron resonant frequency is determined 
from the measured mercury precession frequency z/Hg by 


^n,Hg — 


7n 


7Hg 


^Hg. 


( 3 ) 


The Ramsey lineshape may then be written as 


N = N 


a cos 



( 4 ) 


where Su = z/n,Hg — 1^1 is the difference between the 
first-order resonant-frequency estimate z^n,Hg and the fre¬ 
quency z^i of the applied r.f. pulses. The linewidth Az^ is 
the width at half height of the central fringe, and is given 
by 


2T + 8y7r~2T’ 

where T is the free precession time and t is the duration 
of each of the 7 r /2 spin-flip pulses. The polarization a is 
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related to the fringe maximum and minimum A^max, A^min 
as 


A^max A^min 
A^max H“ A^min 


(6) 


The phase 


(j) = 


^n,Hg — ^0 
- 

Az/ 


( 7 ) 


incorporates the difference between the true resonant 
frequency z/q the first-order estimate z/n,Hg derived 
from the mercury. This difference arises for several rea¬ 
sons. First, there is the inherent uncertainty in the ra¬ 
tio 7n/7Hg as discussed above. Second, and as also dis¬ 
cussed above, (j) alters because the neutrons were of such 
low energy that they populated preferentially the lower 
portion of the storage volume, whereas the mercury sam¬ 
pled the entire volume uniformly; there was a difference 
in height of the centers of gravity of the two systems of 
a few mm. In the presence of a vertical gradient dBzjdz 
in the magnetic field, the average field sampled by the 
mercury would have been slightly different to that sam¬ 
pled by the neutrons. The normalized frequency ratio is 
therefore related to this height difference Ah by 


B! = ±Ah 


dBzjdz 


(8) 


where the + sign corresponds to Bq downwards. 

Third, as will be discussed below, gradients in hori¬ 
zontal components of the magnetic field can increase the 
neutron precession frequency while leaving the mercury 
unaffected. The rotation of the Earth also has a part to 
play, as do minor effects such as light-induced frequency 
shifts in the mercury. 

These factors, which remained essentially constant 
throughout each run, caused the departure of Ra (Eq. 

from unity. There are, in addition, two factors that 
could provide extremely small cycle-to-cycle variations 
in (j). These are the small mismatch in temporal overlap 
of the mercury and neutron frequency measurements, as 
discussed in Section IIIE[ and a genuine neutron EDM 
or some effect that mimics it. 

The HV was applied in a simple pattern of alternat¬ 
ing sign, with an additional four cycles at E = 0 between 
each 16 cycles of HV at a given polarity. The inital polar¬ 
ity for each run was chosen randomly. In principle, one 
could treat each E value within a run as a distinct data 
set, with its own independent Ramsey-curve fit. Our ap¬ 
proach instead has been to fit a single Ramsey curve to 
each entire run, with the frequency shifts (which are very 
small perturbations, well within the noise) then being 
analyzed on a cycle-by-cycle basis. This has the benefit 
of allowing analysis of runs containing insufficient data 
within each E group for a reliable Ramsey fit. 

While the polarity of the E field could be switched au¬ 
tomatically under computer control, changing the direc¬ 
tion of the magnetic field was a much more onerous pro¬ 
cedure. This was normally undertaken every few weeks. 


Eor each run, typically consisting of several hundred 
measurement cycles, the neutron counts were fitted to 
Eq. Q, as shown in Eig. to provide values of ag 
and an average value of the phase (pg for each of the two 
spin states s. (Eor maximum sensitivity, measurements 
were taken repeatedly at four working points, close to 
the half-height of the central Ramsey valley.) Eor the 
first iteration of each of these two fits, the uncertainty ay 
allocated to each data point was simply the square root 
of the number of neutrons counted in that particular spin 
state. The uncertainty in the fitted mercury frequency, 
which resulted in an uncertainty ax in the x coordinate 
of each data point, was incorporated by calculating an 
“indirect” error aj from the slope of each first-iteration 
fitted curve. 




dx^ 


(9) 


and adding it in quadrature to the original error bar to 
obtain a new uncertainty 


O-tot = (10) 

for use in a second iteration of the fit. This correction 
was, however, generally negligible, as the mercury usually 
measured the magnetic field with much higher precision 
than did the neutrons. 



Phase (rads) 

FIG. 2: (Color online) From [3]. Spin-up and spin-down neu¬ 
tron counts for a single run fitted to the Ramsey curve (Eq. 

g. 


An overall average (weighted as the fit uncertainties) 
0 of the two phases was calculated, giving a common 
phase offset for the entire run and for both spin states. 
Eq. was then inverted to yield an individual phase 
shift S(l)g^i for each data point i and for each spin state: 




(f) — arccos 


- Ns 
Ngag 


( 11 ) 


A phase shift was calculated for each data point by 
averaging the over the two spin states. This aver¬ 
aging helped to remove non-statistical fluctuations in the 














4 


total neutron flux, which could, exceptionally, vary by a 
percent or two during a run. 

Any shift in the neutron resonant frequency from 
the value z/q predicted by the fitted Ramsey curve should 
be proportional to this phase shift 5(j)i\ 

Az/ 

6iyoi = 6(pi -. ( 12 ) 

TT 

In the presence of an EDM, this frequency shift should 
be directly proportional to the strength of the electric 
field. A straightforward linear least-squares fit of the 
frequency shifts SuQi as a function of the applied electric 
field E yielded a value for the measured EDM dmeas, with 
its associated uncertainty, for each run. 

We have recently carried out detailed modeling of our 
storage-chamber geometry using the Opera m finite- 
element analysis package, and we find that the average 
electric field within the volume is 1.1% lower than the 
nominal value determined from the applied voltage di¬ 
vided by the separation between the electrodes. . 

Eor completeness, we note that the statistical uncer¬ 
tainty to be expected from an EDM measurement based 
upon the Ramsey technique is [3] 

(13) 

2aETy/N 

Care was taken to verify that the analysis delivered 
the correct sign of EDM: in particular, with E and 
B fields parallel, an increase in neutron precession fre¬ 
quency characterizes a negative EDM. 


III. APPLIED CUTS 

In this Section we list all cuts applied to the data. No 
systematic dependence of the EDM signal upon any of 
these cuts was observed. After all of these following cuts 
were applied, 545 runs containing 175,217 measurement 
cycles and 2.5 x 10^ neutrons remained. 


A. Manual cuts 

After inspection, 10 runs were rejected for various rea¬ 
sons relating to the failure of hardware components such 
as amplifiers, valves, the HV supply and so on. Some in¬ 
dividual cycles within each run were also cut for similar 
reasons: the majority of these were due to problems with 
the delivery of neutrons. A further 26 runs were rejected 
because they had 19 or fewer measurement cycles, which 
is not sufficient for the HV polarity reversal required for 
an EDM measurement. 

These manual cuts removed about 4% of the avail¬ 
able data, after which the above-mentioned 545 runs re¬ 
mained. 


B. Sub-runs 

During part of the data taking, an adjacent experiment 
was using a superconducting magnet. When turned on 
or off, Ra changed by typically 2-3 ppm. The eleven mea¬ 
surement runs during which this occurred were therefore 
split into sub-runs for separate analysis. 

C. Mercury cut 

The x^/z^ distribution from the cycle-by-cycle online 
Hg frequency fitting is shown in Eig. which is repro¬ 
duced from Eig. 10 of [3]. The red “expected” distribu¬ 
tion curve is based upon the model of a perfectly constant 
frequency, whereas in fact we do expect some slow vari¬ 
ation over time within a cycle. That being the case, a 
large jv is not necessarily indicative of a problem with 
the frequency measurement, and we therefore retain a 
generous portion of this distribution beyond the vicinity 
of the peak. A cut was made at jv = 3, by which time 
the tail is fairly flat. 

This cut removed nearly 7% of the remaining data. 

We note in passing the discontinuity at x^/^ = 4. This 
arises because, beyond this point, the online fitting pro¬ 
cedure attempted to correct for potential hardware errors 
(e.g. gain saturation or a missed ADC reading), as dis¬ 
cussed in [3]. 



0 1 2 3 4 5 

Hg xh 


FIG. 3: (Color online) From [3]. The distribution of y^ jv 
from mercury frequency fits. See text for details. 


D. Mercury frequency uncertainty 

The fitted Hg frequency sometimes had a large uncer¬ 
tainty, particularly if the depolarization time was short. 
The distribution of these uncertainties is shown in Fig. 13 
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of [3], and again here on a semi-log plot in Fig. El A typ¬ 
ical value is 1-2/iHz, which (when scaled with 7n/7Hg) 
corresponds to approximately a factor of 5 better than 
the typical inherent neutron frequency uncertainty from 
counting statistics. A cut was made at 25 /iHz, at which 
point the relationship is approximately inverted, with the 
mercury uncertainty entirely dominating the frequency- 
ratio measurement. Beyond this, the data are so impre¬ 
cise that they make essentially no contribution at all. 

This cut removed 4% of the remaining data. 
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the 15 s window, i.e. about 1/20 of the way through the 
Ramsey measurement. The mercury frequency-jump dis¬ 
tribution was truncated at ±60 /iHz, corresponding to a 
change in neutron frequency of 230/iHz over this 1/20 
of the Ramsey period. When averaged over the entire 
Ramsey period, the neutrons would therefore see a fre¬ 
quency shift of up to 11 /iHz, or 0.4 ppm, that would not 
be compensated properly by the mercury. Including the 
aforementioned 1/7 probability, this corresponds (even 
at the extremes of the frequency-jump distribution) to a 
potential error in R' of 0.06 ppm on the rare cycles within 
which such jumps occur, to be compared with a typical 
statistical uncertainty on the neutron frequency of about 
0.7 ppm. 

The frequency-jump cut removed 3% of the remaining 
data. 


F. First-cycle cut 

The first cycle of any run is different from any of the 
others, as the neutron trap and guides are initially empty; 
for other cycles there is likely to be some remnant pop¬ 
ulation from the previous cycle. In consequence the first 
cycle often has an anomalously low total neutron count. 

This cut removed 0.2% of the remaining data. 


FIG. 4: (Color online) The distribution of uncertainties from 
mercury frequency fits. See text for details. 


E. Magnetic-field jumps 

The distribution of Hg frequency jumps, i.e. the dif¬ 
ference in Hg frequency between a given cycle and the 
previous cycle, is shown in Fig. 14 of [3]. There are broad 
tails due to occasional sudden changes in field, for exam¬ 
ple due to the movement of an overhead crane or to a 
mechanical disturbance to the /i-metal shields. 

The mercury and the neutron frequency measurements 
do not have perfect temporal overlap. As discussed in [3], 
the Hg frequency was determined by fitting a 15-second 
averaging period at either end of the cycle in order to 
determine the phase, and hence the integrated phase dif¬ 
ference accumulated over the time between them. The 
corresponding phases for the neutrons, on the other hand, 
are determined by the 2-second spin-flip pulses at either 
end. The field jumps can occur at any point during the 
220 s measurement cycle. If they occur outside the Ram¬ 
sey sequence, they are of little concern. If they occur 
between the 15 s windows, they can be regarded as appro¬ 
priately compensated. That leaves a potential 30 s period 
during which there is a risk of incomplete compensation, 
i.e. 1/7 of such jumps can be expected to affect the mea¬ 
surement to a greater or lesser extent. On average, these 
potentially risky field jumps would occur halfway through 


G. Ramsey residuals outlier cut 

The data points of Fig. reproduced from Fig. 5 of 
[3], show the distribution of stretch values ri of the fits 
to the Ramsey curve: 


where vi is the calculated frequency of the ith cycle, 
is its uncertainty and is the expected frequency for 
that cycle as defined by the Hg frequency, the applied r.f. 
and the Ramsey curve function. Ideally, and in the ab¬ 
sence of any EDM-like signals, this distribution would be 
expected to be a Gaussian of unit width. The continuous 
line is a Gaussian of width 1.06. The true distribution 
departs from this Gaussian at about 4cr. The few points 
lying outside this range tend to be associated with runs 
that have other known problems, for example with inter¬ 
mittent failure of the neutron delivery system. A cut was 
therefore made at ±4cr. Because of the symmetric way 
in which the data are taken, cutting the tails from this 
distribution cannot of itself induce a false EDM signal. 

This cut removed a further 0.05% of the remaining 
data. 


H. High voltage 

A single measurement cycle with a leakage current in 
excess of 60 /iA was removed from the data set. All re- 
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FIG. 5: (Color online) From [3]. Distribution of stretch values 
from the fits to the Ramsey curve 


FIG. 6: (Color online) From [3]. Distribution of the cycle-by- 
cycle average leakage currents I 


maining measurement cycles had leakage currents below 
10 /iA, with the great majority being of the order of a few 
nA; the distribution for both polarities is shown in Fig. 

which reproduces Fig. 20 of [3]. No further cuts relat¬ 
ing to the HV were made: the mercury magnetometer is 
relied upon to compensate for any residual effects from 
this source. Under the assumption that leakage currents 
would run along preferred established paths - rather than 
averaging out over different paths - then if such leakage 
currents were to generate EDM signals, measurement cy¬ 
cles with high leakage currents would be expected to show 
a greater-than-average departure from the Ramsey curve, 
and thus to have high Ramsey-residual stretch values 
(Eg. [l4|). In order to quantify this, the data were binned 
by leakage current, and the distributions of stretch val¬ 
ues within each bin were fitted to Gaussians. The widths 
of these Gaussians are plotted against leakage current 
in Eig. No consistent trend is visible. Eurthermore, 
after processing, no dependence of the measured EDM 
on leakage current was observed, as discussed in Section 
IVIIIEl 



FIG. 7: (Golor online) Departure from the Ramsey curve, as 
measured by the width of the stretch value distribution (as 
in Fig.|^, as a function of the leakage current I. The widths 
are constant to within a few percent, and there is no overall 
trend. 


I. Frequency ratio 


fitted lines. We define the stretch values 


Measurements were undertaken with a range of differ¬ 
ent applied 5-field gradients dB^jdz by pre-adjusting 
currents in field-trimming coils. (Initially, this arose 
through trial-and-error as the system was optimized; we 
then settled upon a more consistent configuration for 
each Bq direction, with a few runs later on used to explore 
the effects of deliberately large gradients.) As discussed 
in Section |IV| these field gradients induced false EDM 
signals. Any nonlinearity in this effect would appear as 
a systematic departure of the points in Eig. from their 


- df 


Zi = 


(15) 


where dmeas^ A is the measured EDM for run i, and 
df. is the corresponding value from the fitted line. The 
distribution of these stretch values is shown in Eig. 
There being no evidence of any correlation, and none 
expected, it was decided not to impose any restriction 
upon the range of Ra. 
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FIG. 8: (Color online) Measured EDM stretch value as a 
function of the frequency ratio Ra 


IV. FALSE EDM SIGNALS 
A. Introduction 

A false EDM signal ^in,Hg,f arises when the trapped 
particles experience a gradient dB^jdz in the presence 
of the electric field E [4]. The effect is easiest to under¬ 
stand in the case of an azimuthally-symmetric field with a 
vertical gradient dBzjdz^ i.e., a slightly trumpet-shaped 
field, as shown in Fig. 


from the side walls of the cell, giving a radial field of 
strength that is proportional to the radius r: 


Br 


r dBz 
2 dz 


(16) 


Consider now a particle moving at speed v that crosses 
the storage cell close to its diameter, as shown in Fig. 


10 As it travels through the electric field, the particle 






FIG. 9: From [4]. Showing the shape of the Bo held lines, 
when there is a positive gradient dBzjdz^ shown in relation 
to an outline of the trap used to store ^^^Hg atoms and ul¬ 
tracold neutrons (UCN). Additional fields having lines that 
both enter and leave through the sidewalls, like the one on 
the right-hand side, do not affect the false EDM signals that 
are generated. 


Because V • 5 = 0, field lines either enter or emerge 


EIG. 10: Erom |3]. Showing the B^y fields (in black) seen by 
a particle going back and forth close to the y-axis. Going to¬ 
wards positive y, the Bxy field rotates steadily anticlockwise 
by about 70° as drawn. The first reflection of the particle 
towards negative y causes an instantaneous anticlockwise ro¬ 
tation by about 110° as drawn. The same two rotations occur 
on the path to, and at, the second reflection. The size of the 
rotations depends on the size of BorlBy. 


experiences in its own rest frame an additional magnetic 
field [12] 


Bv 


V X E 


(17) 


above and beyond the laboratory magnetic field Bq. At 
the start of its trajectory, just left of centre at the bot¬ 
tom, it is subject to the radial field By as well as to the 
sideways By component, yielding a diagonal resultant. 
As it traverses the trap, the By component shrinks and 
then reverses direction, causing a smooth rotation of the 
net additional effective B field. Eventually, the particle 
reaches the far end of the trap, is reflected from the wall 
and begins its trajectory back. However, the By com¬ 
ponent then faces in the opposite direction: therefore, 
after a discrete jump at the point of reflection, the ad¬ 
ditional net effective field component continues to rotate 


























in the same direction. The particle thus sees a rotating 
field in the x-y plane, which, through the Ramsey-Bloch- 
Siegert mechanism, ‘pulls’ its resonant frequency 

away from the central value. When averaged over both 
directions of a given trajectory, the net frequency shift 
is proportional to and it therefore mimics an EDM 
signal. 

The size of this false-EDM effect depends upon the 
relative magnitudes of the orbital-trajectory frequency 
of the particle and its Larmor precession frequency. The 
neutrons are in the adiabatic limit, where their rapidly 
precessing spins can follow the variations in field. In that 
case, the false EDM is given by [4] 


Su = 


^xy dBz 
SttBqC^ dz 


(18) 


where Vxy is the average transverse particle velocity. The 
origin of this shift is a geometric (or Berry’s) phase [15], 
and therefore - as pointed out explicitly in [16] - is inde¬ 
pendent of the coupling to the magnetic field; the gyro- 
magnetic ratio is absent from this equation. 

The mercury lies in the non-adiabatic regime, and its 
frequency shift in such a field gradient is [4] 


6iy = 


7^1 dBz^ 

IGttc^ dz 


(19) 


where is the radius of the storage bottle. This has 
been independently measured m- A more general ex¬ 
pression, valid for arbitrary magnetic fields in the non- 
adiabatic limit, was derived more recently m- 


Si' = {xB^ + yBy) E, (20) 


where the brackets refer to the volume average over the 
(arbitrarily shaped) trap. 

Since the mercury is used to compensate for shifts in 
the magnetic field, any EDM-like component contribut¬ 
ing to the mercury frequency will affect the measurement 
of the neutron EDM. We note that contains a sign 
error in its Eq. 20 and Eig. 2, arising from this transferral 
of the mercury false EDM to the neutron signal. 

In this experiment, the contribution to dmeas of the 
false-EDM effect in the mercury is about 50 times larger 
than the geometric-phase induced false-EDM effect of the 
UCNs. 


B. False-EDM analysis: First iteration 

Combining Eq.j^with Eq.[^ the mercury’s false-EDM 
contribution to dmeas is seen to be 

C^n.Hg.f = ±^l7n7ggl ^^°2 ' R' = ■ R', (21) 

where is the trap radius and the + sign again corre¬ 
sponds to Bo downwards. The notation dn,Hg,f is drawn 
from Eq. 87 of [4]. It follows that we can write 

dmeas = d^ + <^n,Hg,f = d'^ ± k ’ {R' - Rq), (22) 


where d'^ is the true dn plus all other systematic ef¬ 
fects discussed below, and Rq is the value of R' where 
dBz/dz = 0. Eq. 22 defines two straight lines, one with 


a positive slope for Bq down and one with a negative 
slope for Bq up. Naively, one would expect that the 
crossing point (i?x, dx) would occur at dB^/dz = 0, and 
would therefore provide an estimator of d!^ free of dn,Hg,f. 
In reality, as we discuss below, various effects can induce 
shifts in this crossing point, and we therefore must apply 
appropriate corrections. 

As discussed above, Eig. shows the data (binned for 
clarity) for dmeas as a function of R! for each direction 
of Bq. The solid red lines represent a least-squares fit 
to all 545 of the (unbinned) run results, using as free pa¬ 
rameters the two intercepts c \, C 2 and a common absolute 
slope k. This was done by minimizing 


x" = E 


i i kR' Cj 


(23) 


where the terms in the sum that correspond to Bq down 
(up) used the + (-) sign and the intercept ci ( 02 ). The 
fit yielded a crossing point at R^ = 3.50 ±0.81, dx = 
—0.66 ± 1.53 X 10“^^ ecm, with ^ jv = 651/542 and 
k = (1.88 ±0.25) X 10“^^ e cm/ppm. This fitted value of 
the slope was particularly influenced by one run (number 
1900) that was taken at a large applied dB^/dz (giving it 
significant leverage) and which departed from the fitted 
line by 2.7a: excluding that run would have changed the 
value of the fitted slope to 1.63 ±0.27 x 10“^^ e cm/ppm 
and increased the crossing-point EDM value by 0.22 x 
10“^^ e cm, but there appeared to be no a priori reason 
to do so. The run was therefore included both in the 
initial 2006 analysis [T] and in the analysis that follows 
here. 

According to Eq. the slope k is expected to be 
inversely proportional to Ah. We shall discuss the esti¬ 
mation of Ah in Section IV C[ The relationship between 
R and dn,Hg,f is also affected substantially by the phe¬ 
nomenon of gravitationally enhanced depolarization m, 
and this will be discussed in Section IIV PI below. In ad¬ 
dition, the slope k can be altered by a few percent (al¬ 
though still remaining highly symmetric under Bq rever¬ 
sal) by various other mechanisms including the UCNs’ 
own false-EDM signal (a 2% effect); a slight reduction in 
mean free path due to cavities and grooves in the elec¬ 
trodes as well as to the presence of 10“^ torr of He gas 
to prevent sparks [4l[T9] (1%); and a possible bias in the 
volume-averaged frequency measurement due to asym¬ 
metric surface relaxation of the Hg, e.g. if there is a pref¬ 
erential depolarization on one or other electrode (up to 
5%). 


C. Estimating Ah 

As seen in Eq. [^ the slope k of the lines in Eig. 
is determined, to first order, by the height difference Ah 
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between the centers of mass of the mercury and neutron 
ensembles. Initial estimates[T] yielded Ah = 2.8 ±0.1 
mm. This was based on a misinterpretation of the fre¬ 
quency response in a trap with a variable-height lid, aris¬ 
ing from a lack of understanding at the time of the way 
in which gravitationally enhanced depolarization affects 
the measured neutron-to-mercury frequency ratio. A de¬ 
tailed Monte Carlo simulation has now been carried out, 
taking into account the known properties of the neutron 
source, guide tubes, elbow, polarizer foil, and storage 
chamber; it reveals that the spectrum was considerably 
softer, i.e. of lower mean energy, than had originally been 
assumed. 

The simulated spectrum is also essentially consistent 
with a largely analytic calculation [3], which concluded 
that at the start of the storage period the UCN density 
pe{z) at the level z = 0 of the lower electrode could be 
roughly approximated by a softened Maxwell spectrum, 

Pe(0) =A-e^ exp {-e/Vf ), (24) 

where Vf = 88 neV is the Fermi potential of the quartz 
insulator, above which energy p drops to zero. Hence¬ 
forward, when we refer to the UCN energy, it is to be 
taken to mean the kinetic energy of the UCN at the 
level of the bottom electrode; or, equivalently, the to¬ 
tal (kinetic plus potential) energy at any point in the 
bottle, where the potential energy is again referenced to 
the level of the bottom electrode. Eq. estimates rela¬ 
tive densities at the bottom electrode, rather than total 
numbers in the bottle: the latter would have a slightly 
different distribution, since higher-energy neutrons ex¬ 
tend further in height. In order to account for this, we 
note that for a mono-energetic group having energies be¬ 
tween e and e ± de, where the energy e is given in terms 
of the maximum height that the neutrons could reach if 
travelling vertically upwards, the available phase space 
dictates that the density variation with height pe{z) has 
the form [20j 


Pe{z) 



(25) 


and one must integrate this function over the bottle 
height (or, if less, the attainable height e) to establish the 
relative numbers of UCN of each energy trapped within 
the bottle. 

Although our best estimate of the spectrum is that 
given by the simulation - which will underlie the analysis 
that follows - we have carried out a complete analysis 
for both spectra, in order to see whether there was any 
significant sensitivity of the final result upon the initial 
spectrum. These two initial spectra are shown as the 
uppermost line and set of points in Fig. EH 

The spectrum softens during storage, since higher- 
energy UCN not only have a higher interaction rate with 
the walls; they also have a higher loss probability per 



Height-equivalent energy (m) 


FIG. 11: (Color online) Calculated energy spectrum of the 
stored UCN, both at the start (t = 0) of the storage pe¬ 
riod (blue dotted line: softened Maxwell spectrum; purple 
diamonds: simulation) and after t = 135 s (solid red line: 
softened-Maxwell spectrum; blue circles: simulation). The 
total area under the simulated spectrum at t = 0 s is nor¬ 
malised to one; the softened-Maxwell spectra are scaled to 
match the simulated peak height at t = 135 s. 


bounce, as given by m 


m) = V 





(26) 


where f = W/V is the ratio of the imaginary (VF) to the 
real (U) part of the material’s Fermi potential. In fact, 
although V is in general well known, W is difficult to 
determine; and in any case, losses are likely to be domi¬ 
nated by hydrogen that has diffused into the containing 
surfaces. This hydrogen can - because of its extremely 
high incoherent scattering cross section - substantially 
influence loss rates without significantly altering the sur¬ 
face potential V. The value of / used in the simulation 
was therefore adjusted until, at / = 4x 10“^, the loss rate 
(including (3 decay) more or less matched the observed 
loss rate of neutrons in the bottle, as shown in Fig. 
This is broadly similar to the value / = 3 x 10“^ ob¬ 
tained from an independent data-simulation comparison 
undertaken with the apparatus circa 2006 [22]. It was 
assumed in these calculations (as in [22]) that / would 
be similar for both the insulator and the electrodes. The 
softened-Maxwell spectrum was propagated in a similar 
manner (although for this / = 3 x 10“^ gave a better fit 
to the storage-time data, as shown also in Fig. [^. The 
resulting two spectra after 135 s of storage, by which time 
the Ramsey sequence is complete, are shown along with 
the initial spectra in Fig. [^ The resulting final simu¬ 
lated spectrum is slightly firmer than the final softened- 
Maxwell spectrum. 

Throughout this analysis the storage trap was modeled 
as a simple cylinder, with an additional cavity at the 
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bottom centre, 4.0 cm deep and 3.4 cm in radius, within 
which is set the UCN entrance door to the trap. 



FIG. 12: (Color online) Fraction of neutrons remaining in the 
bottle as a function of the time for which they are stored. Red 
triangles represent measurements; blue circles and the solid 
green curve represent retained fractions for the simulated and 
the softened-Maxwell spectra, respectively. 


Given that the calculated Ah after storage is 3.7 mm, 
to first^proximation the anticipated slope k of the lines 
of Fig. [^should decrease from 1.57 x 10“^^ ecm/ppm (in 
m) to 1.18 X 10 ecm/ppm, before taking gravitation¬ 
ally enhanced depolarization into account. 

D. Gravitationally enhanced depolarization 

Under the influence of gravity, neutrons of different 
energies effectively sample different regions of the storage 
trap. The (sometimes surprising) consequences of this 
effect have only recently been studied in detail [iHlEai, 
and have now been validated by comparison with data 
[24l |25] . This stratification results in a dephasing of the 
ensemble in addition to that arising from the intrinsic 
depolarization that naturally occurs within each energy 
bin. Since lower-energy neutrons preferentially populate 
the bottom of the bottle, this dephasing is asymmetric, 
and results in a (strongly spectrum dependent) frequency 
shift. Calculations in m were based upon a Maxwellian 
velocity distribution with a sharp 93 cm height-equivalent 
energy cutoff. The simulated spectrum shown in Fig. 
is much softer than this, showing a clear peak at 20- 
30 cm height-equivalent energy. Under these conditions 
the gravitational-depolarization effect distorts the shape 
of the crossing lines, significantly enhancing the slope 
for measurements beyond a few ppm from the origin, as 
shown explicitly in Fig. The solid red line is the 
best linear fit to the data, mirroring the lines shown in 
Fig-i whereas the dashed green line is the expectation 
from Eq. The dotted blue line includes the effect of 


gravitationally enhanced depolarization; we designate its 
functional form as 

C^n,Hg,f = (27) 



FIG. 13: (Golor online) Plot showing the anticipated false 
edm dn,Hg,f ^ function of R' , for Bo upwards. The dashed 
green line is the expectation from Eq. |21| and the blue dotted 
line shows the revised expectation when gravitational depo¬ 
larization is taken into account. The red line is the linear best 
fit to the data, as shown in Fig. 

The data of Fig. [^were fitted to 

corneas = ± - i?x), (28) 

thus allowing the crossing-point coordinates Rx^dx to 
vary, as well as allowing ^ to be multiplied by a factor 
k' as a consistency check of the input spectrum, since 
the slope of the function is completely determined by the 
spectrum via Ah. As well as crossing-point values of 
Rx = 2.62±0.96ppm, dx = (-0.55T1.51) x ecm, 

the fit yielded k' = 1.15±0.16; and, when it was repeated 
without inclusion of the aforementioned possibly anoma¬ 
lous datum of Run 1900, we found k' = 0.93±0.19. These 
numbers are consistent with unity, giving confidence in 
the simulated spectrum (for comparison, the softened- 
Maxwell spectrum yielded k^ = 1.37 ± 0.19), and so, go¬ 
ing forward, we set k' = 1. We therefore fitted the data 
instead to 

dmeas = ± “ Rx)‘ (29) 

This yielded crossing-point values of Rx = 2.11± 
0.89ppm and dx = (—0.59 ± 1.53) x 10“^^ ecm. jv 
had a value of 650/542 = 1.19, unchanged from the ear¬ 
lier fit. We note that the change from three to two free 
parameters was accommodated almost entirely by shift¬ 
ing Rx , with virtually no change in dx\ we deduce that 
the crossing-point value of the measured EDM dx is rela¬ 
tively insensitive to such changes in the slope, and there¬ 
fore also to the detailed shape of the input spectrum. 
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By eye it may appear that the straightforward linear 
fit to the data is so close to the predicted curve that it 
should make little difference to any analysis. However, 
the central linear region of the predicted curve is now 
shallower by a factor 1.18/1.88 = 0.63 compared to the 
slope of the original linear fit, and any systematic shifts 
in R will therefore have a correspondingly reduced effec¬ 
tiveness in producing systematic shifts in dx • 

A number of assumptions were made while calculating 
the size of the predicted false-EDM effect - for example, 
the gravitational-depolarization calculation here assumes 
only a uniform vertical H-field gradient; and no allowance 
has been made for the effect of grooves in the electrodes 
(which are necessary to locate the insulating trap walls). 
Effects of intrinsic depolarization were explicitly included 
in this model, although perhaps not perfectly since (as 
discussed in [18]) it has a significant dependence upon 
the specularity of reflections within the trap. However, 
its contribution to the frequency-shift effect is negligible 
[24] : primarily, it simply smears the phases. Addition¬ 
ally, it is worth noting that in principle the effective UCN 
spectrum could be slightly influenced, leading to different 
slopes A:^, for the two Bq field directions, if the spin 
relaxation were substantially different for the two config¬ 
urations. This extra relaxation, if it existed, would be 
expected to be UCN-velocity dependent, and so would 
change the effective UCN energy spectrum. However, 
the observed polarisation product a (i.e. the visibility 
of the Ramsey fringes) is identical at the 2% level for 
the two H-field polarities. This suggests m that the 
average UCN energy, and therefore Ah (which is approx¬ 
imately inversely proportional to the energy) should be 
the same to better than a few percent between the two 
polarities. It is difficult to conceive of a mechanism that 
would change the slope more significantly than this upon 
field reversal. 

We have one further small correction to make at this 
point. Since the Ramsey-resonance measurements were 
carried out with an oscillating rather than a rotating 
transverse RE field, the measured frequencies are sub¬ 
ject to a Ramsey-Bloch-Siegert (RBS) shift [131 (Ill ES] 
given by 


E. Statistical sensitivity 


Eor the majority of the data-taking runs, the upper 
electrode was held at voltages of circa ± 80 kV, giving 
an applied electric field of £" = 7 kV/cm. (Voltages as 
high as 130 kV were sustainable towards the end, whereas 
in early runs it was not possible to go much above 60 
kV.) The average polarization was a = 0.58, and the 
Ramsey coherence time was T = 130 s. Removing from 
consideration the cycles measured at E = 0 kV/cm, Eq. 
13 yields an anticipated sensitivity of 1.34 x 10“^^ ecm. 
Our achieved sensitivity is a few percent larger than this, 
for a number of reasons: for example, field changes meant 
that not all cycles were taken at the intended points on 
the Ramsey curve; some cycles had frequency errors en¬ 
larged by relatively rapid depolarization of the mercury; 
and any changes in field gradient resulting in changes in 
R' during the run (which at the ^ 0.2 ppm level or be¬ 
low would not be detectable) could also add noise at the 
percent level. 


V. MODIFICATIONS TO THE FALSE EDM 
SIGNAL 


As indicated above, there are some processes that can 
displace the crossing point, and thus interfere with this 
technique of removal of false-EDM effects - essentially 
any process that changes Ra and/or dn,Hg,f without con¬ 
forming to the ratio between the two given by Eq. 21, 
and where, in addition, the changes differ with the direc¬ 
tion of Bq. These processes may broadly be divided into 
two categories: those that shift frequencies, and hence 
move the lines of Eig.[^ horizontally, and those that move 
the lines vertically by altering the gradient-induced false 
EDM signal relative to the expectation from {dBzjdz). 
In this section we summarize these effects, before describ¬ 


ing in more detail in Section VI the additional diagnostic 
measurements that were carried out first to character¬ 
ize and then to compensate for them. Other systematic 
effects, which contribute to the overall uncertainty but 
that do not introduce a bias to the data, are considered 
in Section [yim 


^Rrbs = 


TT 

4a;2tT (l + ^) ’ 


(30) 


where ujq is the resonant frequency (182 rad/s), t is the 
period (2 s) for which each of the RE pulses is applied, 
and T is the (130 s) period of free precession between the 
RE pulses. This amounts to a mere 0.09 ppm. 

Overall, therefore, we regard the match between the 
measured and expected curves to be in reasonable agree¬ 
ment. Going forward, we use the (RBS-corrected) value 
Rx = 2.02 ± 0.89 ppm, along with dx = (—0.59 ± 1.53) x 
10“^^ ecm, as input to our subsequent analysis. 


A. Uncertainty in 7n/7Hg 

An offset in 7n/7Hg shifts the two lines of Eig. [^side¬ 
ways in the same direction by the same amount, leav¬ 
ing dmeas Unaffected. This therefore has no direct effect 
on the final result, although it can influence slightly the 
overall fit to the data. 


B. Horizontal quadrupole fields 

Eields with finite dB^/dy and/or dBy/dx but with 
{dBxjdx + dBy/dy) = 0 = —dBzjdz cause Ra to in- 
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crease quadratically [4] without contributing to dn,Hg,f. 
This arises because the neutron spins follow the total 
field direction adiabatically, with a correspondingly in¬ 
creased precession frequency, whereas the (nonadiabatic) 
Hg atoms average out the transverse components and 
are generally sensitive only to the Bz component. We 
consider here a quadrupole field aligned with z, with 
Bx = qy, By = qx. This gives rise to a shift in Ra of 

q^r^ 

SRq = (31) 

where, as before, r is the radius of the trap. 

As discussed in [3], the region in the vicinity of the 
EDM measurement volume was scanned with a flux- 
gate magnetometer. This revealed the presence of such 
horizontal field components, quadrupolar in form and a 
fewuT in magnitude throughout the 12 cm height of the 
measurement bottle. The magnitude of these quadrupole 
fields also depends slightly upon the orientation (upwards 
or downwards) of the Bq holding field. This is to be ex¬ 
pected, since the Bq return flux passes through the in¬ 
nermost mu-metal shield: indeed, about 40% of the Bq 
field arises not from the Bq coil directly but from the 
magnetization of the shields. The geometry of the 1.15 
m diameter /i-metal shield around the Bq coil has irreg¬ 
ularities on the scale of a few mm, so it seems reasonable 
to expect that aberrations in the Bq field near the outer 
edges of the trap should be of the order of a few uT. Un¬ 
like the remanent fields arising from any permanent mag¬ 
netization of the shields themselves, any B^xy fields will 
reverse when Bq is reversed. Contributions from these 
two different sources may therefore add in one field di¬ 
rection and subtract in the other, leading to a differential 
shift of all Ra values, and thus of the two lines, thereby 
changing the EDM value dx at the crossing point of Eig. 

m 

Eurther evidence for the existence of such quadrupole 
fields came some years after the EDM measurements de¬ 
scribed in this paper, when the shields were disassem¬ 
bled, shipped to the Paul Scherrer Institut (PSI) and 
rebuilt. The results of a detailed fluxgate scan taken in 
2010 within the rebuilt shields is shown in Eig. and 
here a quadrupole-type field pattern can be seen clearly. 
A detailed harmonic analysis of the field shape yielded an 
average q value of -1.29 (3.24) uT/m for Bq up (down), 
which would imply corresponding quadrupole shifts of 
5Rq = 0.02 (0.14) ppm for the two respective Bq direc¬ 
tions. We do not include these figures in our analysis at 
this point, but we will return to the issue briefly during 
our conclusions. 


C. Rotation of the Earth 

The rotation of the Earth shifts all of the frequency 
ratio measurements Ra to lower values by 1.33 ppm when 
the Bq field is upwards, and to higher values by 1.33 ppm 
when the Bq field is downwards lais]. This therefore 



FIG. 14: (Color online) Results of a detailed scan of the 
magnetic field inside the shields, taken some years after the 
EDM measurements described in this paper. By this time the 
shields had been disassembled, moved to PSI and reassem¬ 
bled. The quadrupole-field pattern is clear, and is obviously 
associated with the shields themselves. The scales on the x 
and y axes represent the horizontal coordinates, in mm. The 
magnitude of the field is represented both by the lengths of 
the arrows and, equivalently, by their colors, as shown in the 
scale (in nT) at the bottom. 


acts in a manner very similar to that of the horizontal- 
quadrupole shifts. 

D. Localized losses 

5-field averaging in the trap is affected by localized 
loss of UCN and Hg particles, and by polarization loss in 
the presence of the 10“^ fractional Bq inhomogeneities, 
which may change with Bq direction. However, we es¬ 
timate that the resulting Ra shifts are < 0.1 ppm and 
< 0.01 ppm for the UCN and Hg respectively. Eurther- 
more, they will be indistinguishable from the quadrupole 
shifts. We do not consider them further. 


E. Magnetic dipoles 

As demonstrated by Eq. the field of a permanent 
magnetic dipole (PMD) close to the trap results in a non- 
uniform dBzjdz^ and therefore induces a frequency-ratio 
shift ^5dip, moving the lines of Eig. horizontally by 
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equal and opposite amounts for the two directions of Bq. 
It also adds to dn,Hg,f a false EDM ddip [27], thus shifting 
the lines vertically by this amount. 

Our fluxgate magnetometer surveys of the trap cannot 
rule out PMD fields of less than 1 nT at 2 cm from the 
inner surface. Large areas of the trap are Si02 or Al, 
backed by large voids, and do not come under suspicion; 
but the mercury and UCN doors involve a heterogeneous 
collection of small parts close to the trap. 

Based upon our simulations [i [22 (since confirmed 
to be perfectly consistent with the analytic prediction of 
(iTj), we included in our 2006 analysis a ddip uncertainty 
of ±6.0 X 10 e cm to allow for an undetected 1 nT 
PMD at the mercury door. 

In the case of the UCN door we have better diagnos¬ 
tics. As discussed above, this door sits at the bottom 
of a small cavity, 4.0 cm deep and 6.8 cm in diameter, 
at the center of the lower electrode. As detailed in Sec¬ 
tion measurements were made of the neutron and 
mercury response in a storage trap with a ceiling of vari¬ 
able height, and the results provided clear evidence for 
the existence of both dipole and quadrupole fields. This 
conclusion was reinforced by measurements undertaken 
some years later when, as part of the process of charac¬ 
terizing the system during its removal to PSI, scans were 
made of the components. A crank-and-sliding-plate as¬ 
sembly that converted rotation of the driving shaft into 
the linear motion of the neutron door was revealed to be 
slightly magnetic, with summary records showing that 
a sensor several cm away registered a field of 2.5 nT. A 
disk that constituted the lower part of the door itself 
also yielded a total field of about 100 pT at a distance 
of ^ 13 cm from the sensor. Various other parts, des¬ 
ignated “screws”, “piston”, “Al frame” and so on, were 
tablulated as having associated fields of a few pT to a 
few tens of pT. Detailed records of orientations and posi¬ 
tions were not kept, since the primary goal at that time 
was merely to determine whether or not the components 
were magnetic. It is therefore not possible at this point 
to specify the dipole position and strength precisely from 
the information now available, so as in the 2006 analy¬ 
sis we rely upon the contemporary measurements and 
regard the later scans merely as qualitative supporting 
evidence, but it seems impossible to avoid the conclusion 
that there were dipole-like fields of ^ nT strength in the 
vicinity of the UCN door. The picture is also consistent 
with measurements undertaken (after the EDM measure¬ 
ments, but prior to moving the apparatus to PSI) with 
a fluxgate magnetometer located in the vertical neutron 
guide immediately underneath the neutron door. This 
registered a slight change (circa 1 nT) in the ambient 
magnetic field when the door mechanism was operated. 

The dipole-field analysis that we have carried out as¬ 
sumes a vertically orientated dipole centered on the axis. 
There is no particuar reason to assume that any PMD 
would have this particular orientation or precise radial 
position. However, we note that the associated frequency 
shifts arise via {dBzjdz)^ which would be much lower 


in the case of a horizontal orientation (not least due to 
cancellation of the respective field components from op¬ 
posite ends of the dipole). Eurthermore, the analysis of 
[TT] (and in particular Eig. 2 therein) shows both that a 
horizontal alignment (for a dipole within a few cm of the 
axis) yields a much smaller EDM contribution than does 
a vertical alignment, and also that the EDM contribution 
(of the vertically aligned dipole, again within a few cm 
of the axis) has very little sensitivity to the radial posi¬ 
tion. We therefore conclude that our model accounts in a 
reasonable manner for the most important contributions. 

The post-hoc magnetic-field scans at PSI revealed no 
magnetic contamination in the vicinity of the mercury 
door, although a re-scan in 2011 showed a small contam¬ 
inant outside the storage volume, at a radius of 31 cm, 
within the threaded hole connecting the ground corona 
ring to the lower electrode. If present throughout our 
measurements, this could have contributed a false EDM 
of up to ±4.0 X 10“^^ ecm (according to the Pignol- 
Roccia analysis). In order to be conservative, though, we 
retain the full ±6.0 x 10“^^ ecm contribution to the sys¬ 
tematic uncertainty, in order to accommodate both the 
possible undetected presence of a (further) PMD as well 
as the fact that we have not modelled the actual door 
dipole perfectly. 


F. Mercury light shift 

The presence of the mercury reading light can, via the 
RBS mechanism nain], shift the resonant frequency of 
the mercury atoms [28l l29] . Such shifts are produced by 
any small component, parallel to Bq, of the ^^^Hg probe 
light beam passing through the precessing ^^^Hg atoms. 
This component, and the consequent Ra shift, reverse 
sign on reversal of Bq. Ealse EDM signals can then arise 
in two distinct ways: 

• Direct: If any changes in intensity are correlated 
with the electric field direction, the resultant fre¬ 
quency shift would mimic an EDM. 

• Indirect: Even if the intensity of the light is com¬ 
pletely constant, the light shifts move the data 
points of the two lines in Eig. in opposite (hor¬ 
izontal) directions, and thus contribute a vertical 
offset to the crossing point. 

A slight dependence of Ra on the incident light inten¬ 
sity was indeed found during the 2006 analysis, the mag¬ 
nitude ^ 0.2 ppm being in agreement with theory. (Sim¬ 
ilarly small shifts were found in a recent measurement 
of the ratio of neutron to ^^^Hg gyromagnetic ratios, us¬ 
ing the upgraded nEDM apparatus at PSI with similar 
Hg discharge bulbs but with a new storage chamber [9].) 
A correction to dmeas was made on a run-by-run basis, 
leading to an overall correction of (3.5 ± 0.8) x 10“^^ 
ecm. Eor this analysis we have simply applied exactly 
the same set of (almost imperceptible) shifts to the data. 
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prior to the crossing-lines fitting procedure. Since the de¬ 
pendence should thereby have been removed, we expect 
no remaining net bias from this source. Details of the 
light-shift analysis are discussed extensively in [3], and 
will not be addressed further here. 


VI. AUXILIARY MEASUREMENTS 

Separate ^Hg and Ra measurements (without E 
fields) were made in an auxiliary trap with a roof that 
could be raised or lowered to change the height H. Criti¬ 
cally, the trap was built on the same lower electrode and 
door mechanism that were used for EDM data taking. 
Full details of the apparatus and of the measurements 
made are available in [30] . 

The trap had a smaller radius than the main data- 
taking bottle, 18.5 cm rather than 23.5 cm, and because 
the floor of the trap consisted of an aluminum plate that 
rested on the lower electrode, the door-cavity depth be¬ 
came 6.0 rather than 4.0 cm. These differences would 
have reduced the quadrupole shifts by a factor Q ^ 0.59. 
(Note that, although the quadrupole shifts are reduced in 
this way, the shifts due to Earth’s rotation are not.) The 
roof and floor were coated with deuterated polystyrene 
(DPS) rather than the diamond-like carbon used for the 
main data-taking trap, but the DPS tended to flake away, 
leaving some bare aluminum exposed: the resulting UCN 
spectrum would therefore have had a slightly lower end¬ 
point cutoff than in the data-taking trap, although after 
some tens of seconds of storage this would make little 
difference. Measurements were taken at storage times 
of several tens of seconds’ duration rather than the full 
135 s, because with the reduced-height bottle the loss 
rates were much higher. Spectrum calculations were car¬ 
ried out exactly as described above for the data-taking 
bottle, assuming the same initial spectrum and per- 
collision loss rates. Appropriate Ramsey-Bloch-Siegert 
frequency-shift corrections [laiii], typically amounting 
to 0.1-0.2 ppm, were applied to each measurement. 

Assuming Bz{z) = bo biz + b 2 z‘^ ^ one can show that 
{dBzldz)v = 0 for a trap height H when z^Hg is the 
same for roof settings at i^/2 and H. This situation was 
approximated by adjusting the trim coils until the same 
z^Hg was obtained for trap heights of 109 and 59 mm. 
With this field established, R' was measured for trap 
heights H of 34, 59, 84 and 109 mm for Bq down. The 
procedure was repeated for Bq up, although in this case 
measurements were only made for trap heights of 59 and 
109 mm. We denote the values of R' in this “gradient- 
compensated” environment as R'^^{H). Naturally, it was 
not possible in reality to trim the gradients perfectly, but 
studying the Hg response enabled appropriate corrections 
to be made. For this analysis we have restricted ourselves 
to the eleven measurement runs for which such gradient 
corrections to R' amount to less than 0.5 ppm. 

The resulting values of R'^^{H) are listed in Table I. 
Their averages for each height and magnetic field con¬ 


figuration are shown in Fig. The original intention 
of this series of measurements had been to determine a 
value for the ratio of gyromagnetic ratios 7n/7Hg- How¬ 
ever, the strong variation with height, and in particular 
the separation of Bq up and down values as the height 
was reduced, led to the conclusion that there was a dipole 
field of strength ^ 1 nT penetrating into the door cavity. 
Such a field could in principle have been produced by 
a pointlike contaminant more or less at the limit of de¬ 
tectability by a fluxgate magnetometer, such as the afore¬ 
mentioned door-actuator crank. It could also in prin¬ 
ciple have arisen from a more dispersed (and therefore 
undetectable) source such as, for example, a slight resid¬ 
ual magnetization within the BeCu door plate or from a 
thermoelectric current as the door mechanism operates 
in vacuo. 
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FIG. 15: (Color online) Values of R' (in ppm) for a “gradient- 
compensated” field (see text) for each Bo-held polarity, as a 
function of the height of the storage volume. Where more 
than one datum was taken for a particular conhguration, the 
point shown represents the average. The curved dotted lines 
represent an approximation to the htted function: they are 
based upon the spectrum calculated for 75 s of storage in a 
bottle of height 84 mm, whereas the ht uses spectra specihc 
to each datum. 


VII. GLOBAL FIT 

In this Section we detail the manner in which each 
contributing component was modeled and incorporated 
into the ht, before moving on to discuss the outcome 
of the ht and the consequent corrections that must be 
applied to the crossing-point EDM measurement 

A. Dipole field 

A simple model incorporating a dipole of moment p 
situated a distance z^[p below the surface of the door was 
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Bo dir. 

H (mm) 

Tmeas (s) 

td ' 

-Ttgc 

(TR 

4 

34 

30 

4.80 

0.73 

4 

34 

50 

3.83 

0.64 

4 

59 

40 

4.73 

1.17 

4 

59 

70 

2.74 

1.03 

4 

59 

70 

2.96 

0.70 

4 

84 

70 

2.45 

0.31 

4 

84 

100 

2.56 

0.22 

4 

109 

70 

2.42 

0.59 

4 

109 

70 

1.60 

0.57 

t 

59 

70 

-2.12 

0.50 

t 

109 

100 

-1.34 

0.18 


TABLE I: Values of i^g^, in ppm, as a function of the direc¬ 
tion of the magnetic field Bo, the height H of the gradient- 
compensated auxiliary trap, and the Ramsey measurement 
time Tmeas- 

found to fit the form of the auxiliary-trap data extremely 
well. The average 2 : component (Bz) of such a field over 
a circle of radius r a distance 2 ; above the dipole is 

^ (4^) (^2 + ,,2)3/2- (32) 

By dividing the trap volume numerically into 1 mm thick 
discs, the additional contribution to Bz from the dipole 
field and from an additional uniform applied gradient 
could be calculated as a function of the height z. Assum- 
ing a spectral distribution of UCN as given in Section 
IV C| (in contrast to the 2006 analysis, which assumed a 
Maxwellian velocity distribution up to the quartz Fermi- 
potential cutoff), the relative height distribution of the 
neutrons in each energy bin was calculated. The contri¬ 
butions to the mercury and UCN frequencies, and hence 
the expected ratio Ra, could then be calculated as a func¬ 
tion of the height H of the trap roof. Each iteration of 
the fit used its selected values for the strength and posi¬ 
tion of the dipole, for and for the quadrupole shifts. 
An appropriate uniform gradient was then applied so as 
to yield the same calculated mercury response at iL = 59 
mm and 109 mm, thus matching the experimental con¬ 
figuration of the auxiliary bottle. The Ra values for each 
bottle height H were then calculated and compared with 
the data before proceeding to the subsequent iteration of 
the fit. We designate the shifts in Ra due to the gradient- 
compensated dipole field as 5R^i^{H). 

B. Quadrupole fields 

The quadrupole-field shifts 6Rq^^ (where the 

subscript arrows indicate the respective Bq directions) 
are assumed to be independent of height (in reasonable 
agreement with the fluxgate measurements), and to be 
proportional to a scale factor Q that depends in princi¬ 
ple only upon the (volume-averaged) square of the trap 
radius, thus yielding Q = 0.59 for the auxiliary trap rela¬ 
tive to the main data-taking trap as discussed in Section 


EH above. 

As noted above, these quadrupole-type fields almost 
certainly arise from the /i-metal shields. It is possible 
in principle that the dependence upon radius may not 
be as expected, particularly towards the outside of the 
data-taking bottle, which extends to about 40% of the 
(58 cm) radius of the inner mu-metal shield. As a test, we 
allowed the Q values for the auxiliary trap to be reduced 
by a further factor /g, which was introduced as a free 
parameter in the fit. Ultimately this made no difference, 
as the fitted value of /g = 0.95 was close enough to unity 
that the final EDM limit was unaltered. We therefore 
omitted this from the final analysis, i.e. we set /g = 1 . 

The quadrupole-field contributions 5Rq to the R' val¬ 
ues within the fit were therefore Q • 5Rq^^^ for the aux¬ 
iliary bottle and 5Rq^^i for the main data-taking bottle. 

C. Earth’s rotation 

As noted above, the rotation of the Earth shifts all of 
the frequency ratio measurements Ra to lower values by 
SRe = 1.33 ppm when the Bq field is upwards, and to 
higher values by the same amount when the Bq field is 
downwards. For either the auxiliary or the data-taking 
trap in isolation this effect would be indistinguishable 
from the differential quadrupole shifts, and in the 2006 
analysis it was not taken into account separately. How¬ 
ever, when extrapolating from one trap to the other a 
correction has to be made since the Earth’s rotation shift 
is independent of the trap radius and therefore of Q. The 
shifts 5Re are therefore built explicitly into the fit in this 
analysis. 

D. Literature 7 ratios 

With no applied gradient and with any dipole and 
quadrupole contributions appropriately accounted for, 
the true value of R' should be recovered. The aver¬ 
age of the two independent literature values of iTn/qHgl, 
namely 3.8424574(30) from [9] and 3.8424560(66) de¬ 
rived from [ 6 ] and [ 7 ], is 3.8424572(27). This was used 
as an input estimator of R^. It is equivalent to R' = 
0.31 ± 0.71 ppm. 

E. Crossing point 

By symmetry, the average quadrupole shift 5Rq = 
+ dRQ^)/2 away from the true value of R^ should 
correspond to the crossing point R^ of the lines of Fig. 

m 

Rx — Rj + ^Rq- (33) 

It is convenient at this point to introduce the parame¬ 
ter 


SRqq = {SRQf - SRqi) 


(34) 
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which represents the quadrupole-field splitting. 

We note that, as shields were opened, closed and de¬ 
magnetized many times throughout the four years of data 
taking, the profile of the field could in principle have 
changed during this time - indeed, the same trim-coil 
settings could sometimes yield variations in Ra values of 
±3 ppm. However, since the overall xt of 1.2 for the fit to 
the lines of Fig. was not excessively large, such changes 
were clearly dominated by alterations in the vertical gra¬ 
dient, which would result in movement along the lines 
rather than displacement of data points from the lines. 
This result is not surprising - the largest holes in the 
magnetic shields were on top and underneath, and this 
would have provided the primary access route for ingress 
of external fields. 

We account for these random variations by increasing 
the uncertainty of the datum for the crossing point Rx 
by a factor = 1.09 in the global fit. 


F. Fitting procedure and outcome 

The rather simple underlying model of a dipole and 
a quadrupole field yields a fairly complex functional 
shape that can only be evaluated numerically. The (five) 
parameters that were allowed to vary were the dipole 
strength and z position, the quadrupole shift SR^^ 
and SRqq . The (thirteen) input data points were (a) the 
literature value of (b) the crossing point Rx from 
the earlier fit to the crossing lines, and (c) the eleven 
auxiliary-bottle measurements R'^^ of R' in the gradient- 
compensated fields at various bottle heights H and stor¬ 
age times ts". 


^gc ~ ^7 + ^s) + SRe- (35) 

The fit used the Levenberg-Marquardt algorithm. The 
resulting best-fit function is represented, together with 
the data, in Fig. For clarity, the data in that plot 
are reduced to six points by averaging over each height 
and Bq up/down configuration, and the plotted function 
was generated with the spectrum appropriate to a 75 s 
storage time in an 8.4 cm high bottle. 

The fit yielded output values of R^ = —0.22T0.64 ppm, 
quadrupole-field shifts of SRq^ = —0.5 ± 1.1 ppm and 
SRq^ = 2.7±1.1 ppm, and a small dipole of moment m = 
—4.2 X 10“^ Am^ situated 1.3 cm below the surface of the 
door - quite consistent with its being a part of the door 
mechanism. Such a dipole would contribute a field 1.1 
nT in magnitude 2 cm above the surface of the door and 
0.3 nT at the level of the electrode surface, and as such 
would be at or below the limit of detectability. In fact, 
the fit was fairly insensitive to the dipole strength and 
position, with different starting values of the parameters 
yielding a range of dipole strengths from 2-8 x 10 ^ Am^ 
correlated with a range of positions from 0.8 to 2.2 cm 
below the door. R^, Rx and SRqq, on the other hand, 
were robust. 


We note in passing that the value of Rx = 0.86 ± 
1.08 ppm that emerges is significantly shifted (approxi¬ 
mately Icr) from the input value of 2.02±0.97 ppm. How¬ 
ever, since dx is essentially independent of such small 
changes in i?x, as discussed above, we do not pursue this 
line of enquiry further - what we are interested in is shifts 
relative to Rx • 

The overall of the global fit was 8.9 for the eight 
degrees of freedom. The quality of the fit is therefore 
rather good, particularly bearing in mind the complex¬ 
ity of the fitting function, and suggests that the model is 
both sound and complete - there is no evidence of con¬ 
tamination from fields of a different form. 


G. Field-shift corrections 

In the absence of the dipole field, we would expect the 
R' values at which the vertical magnetic-field gradient 
within the data-taking bottle is zero to be determined 
by the respective quadrupole-field offsets for each of the 
two Bo directions, together with the shifts due to Earth’s 
rotation: 


^ot,l = ^7 ± ^ (36) 

where the upper (lower) sign corresponds to Bq down 
(up). It would therefore be appropriate under those cir¬ 
cumstances to determine the EDM do at either of these 
points rather than at the crossing point (i?x,dx). By 
symmetry, the values Rq^,R'q^ lie a distance 

\SRe - SRqq/2\ = 0.28 ± 0.19ppm (37) 


either side of the crossing point, i.e. at R' = 1.14, 
0.58 ppm respectively. They therefore share a common 
EDM value do- Using our fitted function from Section 


IV D, we find that this shift in R' away from Rx requires 


us to add Sd = (—0.33 ± 0.14) x 10 e cm to our earlier 
crossing-point value dx = (—0.59 ± 1.67) x 10“^^ ecm 
(where we have expanded the uncertainty on dx by 
as we did for Rx)^ yielding do = (—0.92 ± 1.68) x 10“^^ 
ecm as the EDM value corrected for quadrupole shifts 
and Earth’s rotation. 

Eurther corrections must now be made to accommo¬ 
date the dipole field. The formulation of Pignol and 
Roccia m predicts that the dipole characterized by the 
parameters of our fit would produce, in our storage bottle 
with its door cavity, a false EDM of ddip = 0.44 x 10“^^ 
ecm. The crossing lines of Eig. [^are therefore higher by 
this amount than they would otherwise be, and we must 
compensate by subtracting ddip from all of the measured 
data, and thus, ultimately, from dx • 

In addition, this dipole field applies a volume-averaged 
magnetic-field gradient dB^/dz of 0.26 uT/m, which 
shifts the normalized frequency ratio R! by SR^^ = 
±0.98 ppm (where the subscript “dg” indicates “dipole 
gradient”, and the upper sign once again corresponds to 
Bq down). This shift of the lines in opposite directions 
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moves the crossing point downwards by 1.15 x 10“^^ e cm. 
The net dipole-related correction that must be applied is 
therefore now (0.71 ±0.07) x 10“^^ ecm. As mentioned 
above, the dipole strength and position arising from the 
fit are strongly correlated, so the uncertainty on this to¬ 
tal dipole-related shift was calculated by studying the 
behaviour of as a function of the shift for a variety of 
dipole strengths and positions. 

We also at this point add in quadrature to the uncer¬ 
tainty the contribution of 0.6 x 10“^^ e cm from a possible 
mercury-door PMD, as discussed in Section VE| 

Taking all of these dipole-related components into ac¬ 
count, the net EDM thus far - compensated for the 
false-EDM effects arising from field gradients, dipole 
and quadrupole fields, and Earth’s rotation - is dfec = 
(—0.21 ± 1.79) X 10“^^ ecm, where the subscript “fee” 
stands for field-effect compensated. 


H. Consistency check: Polarization data 

Because UCNs of different energies have different val¬ 
ues of Ah, the surviving UCN polarization within a trap 
decreases in proportion to [dB^ldz)^. One would there¬ 
fore expect that a would be maximized close to the B! 
value (for each Bq direction) at which {dBzIdz) is zero, 
as given by Eq. with an additional correction for the 
dipole shift hi?dg- In fact, our calculations show that, due 
to the dipole field, {dBzjdz)^ is minimized at 0.92 ppm 
below for Bq up, and 0.92 ppm above for Bq 
down. We therefore expect the polarization a to be max¬ 
imized at R' = 0.2 ± 1.1 ppm, 1.5 ± 1.1 ppm for Bq up, 
down respectively. 

During runs prior to January 2000 the gradient was 
varied frequently and over a fairly wide range. We have 
plotted in Eig. the polarization a as a function of 
R' for this subset of runs, which represent 14% of the 
total statistical weight of the EDM data analyzed here. 
Although other factors (e.g. loss of polarization during 
filling) affect the data to a greater or lesser extent, the 
envelope of these data points ought to peak in the region 
where the magnetic-field gradient is minimized. 

In the 2006 analysis, bands of points along the tops of 
these distributions were fitted to Gaussians to find the 
positions of the peaks. The fit was carried out multi¬ 
ple times, each time repeatedly removing points that lay 
more than pa below the curve, where p was varied be¬ 
tween 4.0 and 1.0, such that the final remaining points 
would represent the envelope. In recognition of other 
factors contributing to depolarization, and in order not 
to allow the few data points with particularly small un¬ 
certainties on a to skew the fits unduly, an additional 
uncertainty of 1% was added in quadrature to each er¬ 
ror bar. Points with error bars are those used in the fits 
represented by the solid lines. 

Although these results were used in the global fit in 
the 2006 analysis, we feel upon reflection that the large 
scatter of data points underneath each peak casts doubt 
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FIG. 16: (Color online) Polarization product a as a function 
of R', for (a) Bq up, spin down; (b) Bq up, spin up; (c) Bq 
down, spin down, and (d) Bq down, spin up. Solid lines are 
Gaussian fits to points along the upper envelope of the data 
set (which are shown with error bars). 


upon their validity for this purpose. In addition, the fit¬ 
ted peak position has some dependence both upon the 
number of data points included and upon the functional 
form of the fitted curve - which, as we now understand 
[24] . can resemble a sharp peak rather than a Gaussian. 
We therefore present them here only as a consistency 
check, and note that, although the R' values at which 
the fitted curves peak appear to be arguably just a lit¬ 
tle higher than expected, they are close enough to be 
regarded as being in broad agreement. 


VIII. OTHER SYSTEMATIC ERRORS 

We now consider systematic errors that do not involve 
the field-gradient induced false-EDM effect, and that do 
not bias the results. Eor those that are discussed in detail 
in [3], we summarise the nature of the effect and state 
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their contributions to the overall uncertainty. We begin 
with first-order v x E effects. 


A. V X E effects 


The Lorentz transformation of electric and magnetic 
fields to a moving reference frame is m 


B' = -f{B- 


V X E 


0^(7 + 1) 


v{v ' B), (38) 


where 7 = 1 />/! — jc?. As noted above (Eq. 17), in 
the context of the false-EDM effect, a particle moving 
fairly slowly (7 ~ 1) through an electric field E therefore 
experiences an additional magnetic field 


R, ^ — 


V X E 


(39) 


above and beyond the laboratory magnetic field Bq. 

Such fields are clearly linked to the electric field, and 
the necessity of controlling the consequent spurious ef¬ 
fects provides a strong constraint on the design of the 
experiments. In particular, to generate a systematic error 
in the EDM requires that T, E and Bq have components 
that are mutually perpendicular. It was therefore neces¬ 
sary in this experiment to keep E and Bq closely aligned 
and the velocity v averaged over the storage time as small 
as possible. Only the motion of the center of mass of the 
UCN gas and any net rotation of the gas about the center 
of mass can contribute to this first-order vxE effect. The 
effect has been clearly seen in atomic beam experiments 
[ 3 TI 433 ] and has been a cause for concern in neutron beam 
experiments. 

If the neutrons have an average center-of-mass velocity 
of T] perpendicular to E, and E has a component eE 
perpendicular to Bq, the vxE component of the magnetic 
field will give a systematic error in the EDM of 

\d\ = (40) 

= 1.3^6 X 10“^"^ ecm. (41) 


The only source of indefinitely sustained motion that 
can give rise to a finite average velocity would be a net 
upwards movement of the center of mass due to warming 
of the neutron gas as a result of inelastic collisions with 
the walls. The effects of such collisions, in which the 
neutron energy changes by only a small amount so that 
it remains in the UCN range, have been looked for by 
Richardson m- Eor a greased surface, an upper limit of 
50 feV for the energy change per wall collision was estab¬ 
lished, from which it can be calculated that the maximum 
movement will be no more than I mm during the 130 s 
Ramsey measurement period. If the volume-averaged an¬ 
gles between E, B and v are each as high as 2° = 0.035 
radians, the induced false EDM will be I x 10“^^ ecm. 


(Note that such warming is distinct from the changing 
Ah arising from the preferential loss of faster UCN: the 
latter do not ultimately contribute to the EDM measure¬ 
ment, and so cannot generate a systematic error.) 

If there is a difference in a transverse coordinate x 
between where the neutrons enter the trap and their av¬ 
erage center of mass during the storage time, the net 
average velocity is r] = x/Tg. The effect of any such mo¬ 
tion will be greatly reduced by the fact that it is only 
motion during the storage time between the two oscillat¬ 
ing fields that affects the precession frequency. Since the 
trap is filled for 20 s, which in our apparatus is roughly 
1.3 filling-time constants, and this is followed by a 5 s 
settling period, the stored neutrons have spent an av¬ 
erage of one filling-time constant in the trap before the 
first oscillating field is applied. Collisions with the wall 
are expected to destroy any non-random motion within 
100 bounces, i.e., about 10 s, so that any such motion 
that results from the filling process will mostly have died 
away before the first oscillating field is applied. The po¬ 
sition of the center of mass during storage is expected to 
be very close to the axis of the 67 mm diameter guide 
tube. In the unlikely event that the neutrons fill prefer¬ 
entially at X = I cm towards one side of the guide tube, 
and ignoring any damping of the motion in the trap dur¬ 
ing filling, we would be left with a net motion during the 
130 s Ramsey measurement of up to ^ 5 mm. This ef¬ 
fect would then contribute 2 e x 10 “^^ e cm = 6 x 10 “^^ 
ecm. This dominates the earlier result associated with 
warming, and therefore becomes our total uncertainty 
for first-order translational v x E; to be conservative we 
round the value up to I x 10“^^ e cm. 

Higher-order v x E effects, in particular from the net 
increase in |Bo| as the V X E component is added to it in 
quadrature [34], contribute < 10“^^ ecm. 

In a similar manner to the translational effect, any net 
rotational flow of the UCN in conjunction with a radial 
component of the E field may lead to an induced EDM 
signal. However, any such flow of UCN is expected to 
be attenuated by wall collisions before the first Ramsey 
pulse is applied. This systematic effect requires the elec¬ 
tric field to have a radial component cEq, such as may 
exist near the outside of the trap if the insulator of the 
trap becomes charged. If a net fraction /rot of the neu¬ 
trons have such a flow with a velocity of I m/s, and if 
the flow persists for 10 s after the storage period has be¬ 
gun, and hence 5 s into the Ramsey measurement, the 
neutrons will have moved a net ^ 5 m through this field 
within the Ramsey measurement. The effect on the EDM 
is then 5/rot^ x 10“^"^ ecm. It is reasonable to assume 
that only a very small fraction such as /rot ^ 0.001 could 
persist in orbits through this anomalous field region, and 
the radial field fraction e is unlikely to be larger than 10 % 
of the primary field, giving a net error from this source 
of below 5 X 10“^^ ecm. 

There will also be additional cancellations to the v x E 
effect from the reversal of Bq, but in order to be conser¬ 
vative we do not consider these. 
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B. Uncompensated magnetic field fiuctuations 

There may in principle be residual effects from B field 
fiuctuations, such as hysterisis in the /i-metal shield fol¬ 
lowing disturbances in the stabilised Bq coil current sup¬ 
ply caused by pickup from the high-voltage changes. As 
discussed in [3] , we would expect this to manifest itself as 
a dipole-like field originating at the HV feedthrough. 
Because of the effective height difference Ah between the 
centers of the neutron and mercury systems, the experi¬ 
ment retains a slight vulnerability to such a systematic. 
The shielding factor is expected to be about 

^ = — (42) 

dBd 3 Ah’ ^ ^ 

where r ~ 55 cm is the distance from the source of the 
field to the center of the trap. In our 2006 analysis, and 
in [3], we used Ah = 2.8 mm, and thereby concluded 
that the protection factor was about 66. In light of our 
revised calculations of the spectrum, which yield Ah = 
3.7mm, we now reduce this to a factor of 50. 

Fig. 18 of [3] showed the apparent EDM signals of 
the neutron and mercury channels individually, plotted 
against one another. Fig. here shows the same data, 
but this time binned for clarity. 



FIG. 17: (Color online) After [3] Fig. 18. Apparent neutron 
EDM signals (due to uncompensated random magnetic-field 
fluctuations) as a function of the corresponding apparent mer¬ 
cury EDM signals. Data are binned for clarity. 


The neutrons alone yielded a net uncompensated EDM 
signal of (17 ± 4) X 10 e cm. We therefore expect the 
mercury-magnetometer compensation to shield us from 
this systematic effect to a level of (17 x I0“^^)/50 = 
3.4 X 10“^^ ecm. 


C. Mercury atom EDM 

The measured EDM values dmeas are obtained from 
a linear fit to the ratio z^n/^Hg versus E. In principle 


therefore, dmeas contains a contribution from the intrinsic 
EDM ^Hg of the ^^^Hg atom. The true hng has been 
shown to be (0.49 ± I.29stat A 0-76syst) x 10“^^ ecm [35], 
so the systematic error thereby introduced into dmeas is 
a negligibly small (—2 ± 6) x 10“^^ ecm. 


D. Electric forces 

Another possible source of systematic error arises 
from electrostatic forces, which may move the electrodes 
slightly. Discussion of this issue in [3] concluded that the 
associated systematic uncertainty is 0.4 x 10“^^ ecm. 


E. Leakage currents 


If the leakage current that flows through (or along the 
surface of) the insulator between the electrodes has an 
azimuthal component, a component of the magnetic field 
due to the current would be parallel (or anti-parallel) to 
Bq and would produce a frequency shift that changes sign 
when the polarity of the electric field is reversed, giving 
rise to a systematic error in the EDM. As discussed in [3], 
the false signal that would result is likely to be no larger 
than O.I X 10“^^ ecm. As is clear from Fig. 18, which 
is drawn from |3] , no dependence of frequency shift upon 
the leakage current is apparent in the data. 



FIG. 18: (Golor online) From [3]. Frequency shifts (multi¬ 
plied by the polarities of the electric and magnetic fields) as 
a function of leakage current. 


F. Sparks 

High-voltage breakdown within the apparatus can pro¬ 
duce localised high current densities, and can in principle 
lead to permanent changes in the residual magnetization 
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of the magnetic shield. If these sparks occur preferen¬ 
tially for one direction of the electric field, systematic 
effects could be induced in the precession frequency, and 
hence in the EDM signal. A similar effect could be pro¬ 
duced by hysteresis in the innermost shield following a 
disturbance to the Bq supply. However, the mercury 
would naturally compensate for any such effect, just as 
with any other shifts in the magnetic field. 

As sparks invariably disrupt the mercury frequency 
measurement, cycles that contain them are excluded from 
the analysis, so beyond the residual effects just discussed 
the sparks themselves cannot contribute to any artificial 
EDM signals. We therefore do not associate any addi¬ 
tional systematic uncertainty to this effect. 



G. HV AC ripple 

A “ripple” on the high voltage would generate an oscil¬ 
lating displacement current in the storage chamber and 
thereby an oscillating B field. Through the Ramsey- 
Bloch-Siegert mechanism [13], this could in principle lead 
to changes in precession frequency. A detailed discus¬ 
sion of this effect in [3] concludes that such a mechanism 
would in this case not contribute to the overall uncer¬ 
tainty at a level higher than 1 x 10“^^ ecm. 


FIG. 19: (Color online) Stretch values from the fits of Fig. 
in the order in which the data were taken 


by our analysis. No further uncertainty is assigned to 
any of these effects. 


K. Stability of results 


H. Effect of nonuniform mercury depolarization 

The detremental effect of the high voltage upon the 
mercury depolarization time could result in a false signal 
if (a) the average depolarization time were different for 
the two HV polarities, and (b) the mercury frequency had 
some small dependence upon the depolarization time. As 
discussed in [3| , evidence for such an effect was sought in 
the data, and it was ruled out at the level of 1.2 x 10“^^ 
e cm. 


As a final check for possible unforeseen systematic er¬ 
rors, we consider the stability of the measurements. On 
a run-by-run basis, one can consider the stretch values 
(Eq. from the geometric-phase line fit of Eig. [^as a 
function of time, as shown in Eig. 19 this removes the 
natural variation in the measured EDM due to the chang¬ 
ing Ra values. If there were an offset in one series of data 
due, for example, to a magnetic anomaly that was later 
removed, one would see points typically above the line 
for a series of runs at one period in time and below the 
line at other times. No such trend is visible. 


I. Artifacts of the measurement process 

In this class of systematics there is usually no real 
change in precession frequency — only an apparent 
change due to a malfunction of the measurement process. 
An example might be bursts of false counts in the neu¬ 
tron detection channel derived from a high-voltage spark. 
No evidence of any such effects has been seen, and in any 
case any artificial signals from such effects would tend 
to cancel upon reversal of Bq. We do not associate any 
further uncertainty with this mechanism. 

J. Accuracy of Hg frequency measurements 

A number of mechanisms can affect the frequency mea¬ 
surement of the Hg magnetometer. These are discussed 
extensively in [3|, and are already implicitly dealt with 


L. Systematic error summary 

Table H lists all of the systematic errors that we have 
discussed, and gives their uncertainties. Its last line spec¬ 
ifies the net shift in EDM value, and the total systematic 
uncertainty. Eor the avoidance of doubt, we note that 
the shifts listed represent the offsets generated by the 
biases in question. The compensating corrections that 
have been applied are therefore equal in magnitude and 
opposite in sign to the stated values. 


IX. RESULTS AND CONCLUSION 

We have reanalyzed the data that were used in 2006 
to calculate an upper limit on the magnitude of the neu¬ 
tron electric dipole moment [1]. As in the 2006 analy¬ 
sis, earlier data from the interim result of [36] have been 
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Effect 

Shift 

a 

uug light shift (included in dx) 

(0.35) 

0.08 

xt = 1-2 adjustment 

0 

0.68 

Quadrupole fields & Earth rot’n 

0.33 

0.14 

Dipole field 

-0.71 

0.07 

Hg door PMD 

0.00 

0.60 

v X E translational 

0.000 

0.001 

v X E rotational 

0.00 

0.05 

Second-order v x E 

0.000 

0.000 

Uncompensated B drift 

0.00 

0.34 

Hg atom EDM 

-0.002 

0.006 

Electric forces 

0.00 

0.04 

Leakage currents 

0.00 

0.01 

AC fields 

0.000 

0.001 

Nonuniform Hg depol. 

0.000 

0.001 

Total shift of dx 

-0.38 

0.99 


TABLE II: Summary of systematic errors and their uncer¬ 
tainties, in units of 10“^® ecm. Correction for the mercury 
light shift is already incorporated run by run prior to the 
crossing-lines fit; other corrections are then applied to the 
crossing-point EDM value dx- 


Analysis stage 

EDM 

a 

Crossing point dx 

-0.59 

1.53 

Gradient-corrected do 

-0.92 

1.68 

Dipole-corrected dfec 

-0.21 

1.79 

Final result dn 

-0.21 

1.82 


TABLE III: Summary of the net EDM arising from each stage 
of the analysis, in units of 10“^® ecm. 


excluded. We have here taken into account recent de¬ 
velopments, including an understanding of the process of 
gravitational depolarization and more detailed calcula¬ 
tions of the spectrum of the stored UCN, and we have 
also taken into consideration measurements made in the 
meantime on the apparatus. The picture that emerges 
appears to be perfectly self-consistent. 

Once the crossing-point EDM value dx = (—0.59 ± 
1.53) X 10“^^ ecm is corrected for the systematic biases 
listed in Table II, as shown explicitly in Table III for each 
stage of the analysis, we obtain (by adding in quadrature 
the statistical and systematic uncertainties) a final value 
of dn = (-0.21 ± 1.82) X 10-2^ ecm. The 90% and 95% 
confidence limit (CL) ranges are therefore —3.2 < d^ < 
2.8 X 10“^^ ecm and —3.8 < dn < 3.4 x 10“^^ ecm 
respectively. Taking limits that are symmetric about zero 
yields |dn| < 3.0 x 10“^^ ecm (90% CL) and |dn| < 3.6 x 
10-2^ ecm (95% CL). 

This analysis was also repeated with a number of vari¬ 
ations: for example, with the spectrum derived from the 
softened-Maxwell approximation, which yielded (with 
the same uncertainty) central values of 0.86 and 0.79 
X10“^^ e cm with the slope multiplier k' kept fixed and 
allowed to vary, respectively; and without the apparently 
anomalous Run 1900 (resulting in central values of -0.05 


xl0“^^ ecm with the simulated spectrum, and +0.43 
xl0“^^ ecm with the softened-Maxwell spectrum). Fur¬ 
thermore, we have also looked carefully at the implica¬ 
tions of the fields as defined by the maps produced by the 
2010 scans at PSI. We do not consider it appropriate to 
apply the results directly to our data, not only because 
the field scans were carried out after a lapse of several 
years and following dismantling, shipping and reassem¬ 
bly of the shields, but also because the accuracy of mea¬ 
surements of small transverse field components, typically 
a fraction of a percent of the main vertical field strength, 
is limited by the mechanical precision of the mapper sys¬ 
tem itself. Nonetheless, we considered it useful to carry 
out the analysis with the quadrupole shifts 6Rq fixed at 
the values defined by these more recent field measure¬ 
ments. The resulting value of dn was 0.76 x 10“^^ ecm 
(again with the same uncertainty), although as might be 
expected jv for the global fit deteriorated somewhat 
to 1.48. The range of variation of all of these results is 
of the order of the uncertainty on our systematic errors, 
and all lie comfortably within our confidence-limit range. 

Ultimately, and in conclusion, after a thorough re- 
evaluation of the previously known systematic effects 
and the inclusion of various other newly established ef¬ 
fects, the overall corrections to the 2006 result have been 
found - both individually and collectively - to be rel¬ 
atively small compared with the statistical uncertainty. 
As a direct consequence, our new limit represents a de¬ 
parture of only a few percent from the earlier limit of 
jdnl < 2.9 X 10“^^ ecm (90% CL) [T] despite the very 
real differences between the two analyses. 

At the time of writing, the apparatus employed in this 
experiment is enjoying something of a renaissance. It 
has been very substantially upgraded, including in par¬ 
ticular with the provision of an array of high-precision 
cesium magnetometers m, for a further nEDM mea¬ 
surement [38] at the new UCN source at PSI [39]. Fur¬ 
thermore, an innovative new spin-echo based technique 
[25] allows both the measurement of the UCN spectrum 
within the apparatus as well as an accurate determina¬ 
tion of magnetic-field gradients. It is currently running 
with a greater sensitivity than ever before, and, with 
potential systematic errors now being very tightly con¬ 
strained, new results will be limited only by statistics for 
some time to come. 
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